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The effect of strong quantizing magnetic field on the equation of state of matter at the outer crust 
region of magnetars is studied. The density of such matter is low enough compared to the matter 
density at the inner crust or outer core region. Based on the relativistic version of semi-classical 
Thomas-Fermi-Dirac model in presence of strong quantizing magnetic field a formalism is developed 
to investigate this specific problem. The equation of state of such low density crustal matter is 
obtained by replacing the compressed atoms/ions by Wigner-Seitz cells with nonuniform electron 
density. The results are compared with other possible scenarios. The appearance of Thomas-Fermi 
induced electric charge within each Wigner-Seitz cell is also discussed. 
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1. INTRODUCTION 

The theoretical investigation on the properties of compact stellar objects in presence of strong quantizing magnetic 
field have gotten a new life after the discovery of a large number of magnetars [II, H, [H, [H . These exotic stellar objects 
are believed to be strongly magnetized young neutron stars. Their surface magnetic fields are observed to be > 10 15 G. 
Then it is quite possible that the fields at the core region may go up to 10 18 G (the theoretical estimate may easily 
be obtained from scalar virial theorem). The exact source of such strong magnetic field is of course yet to be known. 
These objects are also supposed to be the possible sources of anomalous X-ray and soft gamma emissions (AXP and 
SGR). Now, if the magnetic field is really so strong, specially at the core region, it must affect most of the important 
physical properties of these stellar objects and also some of the physical processes, e.g., the rates / cross-sections 
of elementary processes. In particular the weak and the electromagnetic decays / reactions taking place at the core 
region will change significantly. 

The strong magnetic field affects the equation of state of dense neutron star matter. As a consequence the gross- 
properties of neutron stars [HjflHIH]) e.g., mass-radius relation, moment of inertia, rotational frequency etc. should 
change significantly. In some recent work we have developed the relativistic version of Landau theory of Fermi liquid 
in presence of strong quantizing magnetic field to obtain the equation of state of dense neutron star matter 0, [l(| ■ We 
have also shown that the nucleons acquire complex mass in the a — lo — p meson exchange type mean field in presence 
of strong magnetic field. It has been noticed that due to the complex nature of neutron and proton energies, there is a 
kind of relaxation or oscillation in the iso-spin space. In the case of compact neutron stars, the phase transition from 
neutron matter to quark matter, which may occur at the core region, is also affected by strong quantizing magnetic 
field. It has been shown that a first order phase transition initiated by the nucleation of quark matter droplets is 
absolutely forbidden if the magnetic field strength ~ 10 15 G at the core region [H|, [l2j]. However, a second order 
phase transition is allowed, provided the magnetic field strength < 10 20 G. This is of course too high to achieve at the 
core region. The study of time evolution of nascent quark matter, produced at the core region through some higher 
order phase transition, shows that in presence of strong magnetic field it is absolutely impossible to achieve chemical 
equilibrium (/3-equilibrium) configuration among the constituents of the quark phase if the magnetic field strength is 
as low as B ~ 10 14 G. 

The elementary processes, in particular, the weak and the electromagnetic decays/reactions taking place at the 
core region of a neutron star are strongly affected by such ultra-strong magnetic fields [HI, 03] • Since the cooling of 
neutron stars are mainly controlled by neutrino/anti-neutrino emissions, the presence of strong quantizing magnetic 
field should affect the thermal history of strongly magnetized neutron stars. Further, the electrical conductivity of 
neutron star matter, which mainly comes from free electron gas within the stars and directly controlls the evolution 
of neutron star magnetic field, should also change significantly [14j . 

In another kind of work, the structural stability of such strongly magnetized rotating objects are studied. It has 
been observed from the detailed general relativistic calculation that there are possibility of some form of geometrical 
deformation of these strongly magnetized objects from their usual spherical shapes [l|| [lfl [ljj • In the extreme case 
such objects may either become black strings or black disks. It is quite possible to have gravity wave emission from 
these rotating magnetically deformed objects (see also [HI for the magnetic deformation of nucleonic bags in neutron 
star matter). 
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In a recent study of the microscopic model of dense neutron star matter, it has been observed that if most of the 
electrons occupy the zeroth Landau level, with spin anti-parallel to the direction of magnetic field, and only a few 
are with spin along the direction of magnetic field with non-zero Landau quantum number, then either such strongly 
magnetized system can not exist or such a strong magnetic field can not be realized at the core region of a neutron 
star [19(. We have observed that 10 22 G is the upper limit of magnetic field which the core of a neutron star can 
sustain. We have also noticed that since the electrical conductivity of the medium becomes extremely low in presence 
of ultra-strong magnetic field, the magnetic field at the core region must therefore decay very rapidly. Hence it may 
be argued that the magnetic field at the core of a magnetar can not be too high, it is only the strong surface field 
which has been observed. 

Similar to the study of quark-hadron deconfinement transition inside neutron star core in presence of strong quan- 
tizing magnetic field, a lot of investigations have also been done on the effect of ultra-st rong magnetic field on chiral 
properties of dense quark matter [Ifl, HU, [H, [H, [H, [H, [H, H3, HI, HE HO, HU (see also [H| for the chiral properties 
of strongly magnetized electron gas). The strong external magnetic field acts like a catalyst to generate mass dynam- 
ically. This is one of the most important effect of ultra-strong magnetic field on the properties of charged elementary 
particles. This effect is of course impossible to verify in the terrestrial laboratory. Only the cosmic laboratory offers us 
the opportunity to test this intrinsic phase transition in the world of elementary particles. Therefore, the presence of 
strong magnetic field in stellar matter causes a lot of interesting and significant changes in the properties of elementary 
particles. Also a lot of new phenomena are observed in such cosmic laboratory in presence of ultra-strong magnetic 
field. 

The exotic process of magnetic photon splitting, i.e., the decay of a photon into two photons (or combination of 
two photons into a single one) in presence of a very strong magnetic field, has recently attracted renewed attention, 
mainly because of the great importance this process may have in the interpretation of the spectra of cosmic 7-ray 
burst sources. The basic formulae for magnetic photon splitting had already been derived in the seventies [33| (see also 
[H, HH, Hll ) • Again, since such strong magnetic field can not be generated, even in the pulse form, in the laboratory, 
it is absolutely impossible to verify this exotic phenomenon of photon splitting. 

In this article we shall develop an exact (within the limitation of Thomas- Fcrmi-Dirac model) formalism of relativis- 
tic version of Thomas-Fermi-Dirac model in presence of strong quantizing magnetic field. From this model calculation, 
we shall obtain the equation of state of relatively low density crustal matter of magnetars, which is mainly dense iron 
crystal. We shall replace the iron ions / atoms by Wigner-Seitz cells with non-uniform electron density within each 
cell. We have organized the paper in the following manner: In the next section, we shall present the basic formalism 
of relativistic Thomas-Fermi-Dirac equation in presence of strong quantizing magnetic field. In section 3, we shall 
evaluate the equation of state of low density crustal matter of strongly magnetized neutron star. The kinetic energy 
part of the electron gas within each cell is obtained in section 4. Section 5 is dealt with the electron-nucleus interaction 
energy. In section 6, we evaluate the electron-electron direct interaction energy part. The exchange interaction part 
is obtained in section 7. The possibility of Thomas-Fermi induced charge within each cell is discussed in section 8 
and finally in section 9, the last section, we shall discuss the results and future prospects of this model. 

It is worth mentioning that the formalism we are presenting in this article is also applicable to strongly magnetized 
white dwarfs. For the sake of completeness, we have compared the results of the present model, with other possible 
physical scenarios, e.g., non-relativistic case with zero and non-zero values of Landau quantum number, with the 
non-relativistic and relativistic field free cases, etc. 

Finally, we would like to mention, that the properties of low density magnetized crustal matter, mainly the elec- 
tromagnetic properties, which also includes the transport properties have been studied thoroughly by Potekhin and 
Potekhin et. al. [37j]- 

2. RELATIVISTIC THOMAS-FERMI-DIRAC MODEL IN QUANTIZING MAGNETIC FIELD 

In the past few decades, a lot of work have been done on both Thomas-Fermi and Thomas-Fermi-Dirac models in 
absence as well as in presence of strong magnetic fields to obtain the equation of state of low density crustal matter 
of compact stellar objects. Unfortunately, the calculations are either based on some crude approximation or are 
incomplete in nature. In those calculations, in presence of strong quantizing magnetic field, the occupancy of only 
the zeroth Landau level by electrons were considered. However, such approximation is valid if the magnetic fields 
are extremely strong. In one of the previous calculations we also made a preliminary study of low density crustal 
matter of neutron stars in presence of strong quantizing magnetic field using the Thomas-Fermi-Dirac approach. The 
model was non-relativistic in nature and occupancy of only the zeroth Landau level was taken into account (38j . In 
some work, more than a decade ago, Shivamogi and Mulser did some relativistic calculation for atoms in ultra-strong 
magnetic field [39]. Again, the calculation is applicable for the system where only the zeroth Landau level is occupied. 
A similar type of calculation was also done without magnetic field by Ruffini [401 ] . At this point we should mention 
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that Thomas-Fermi model is not an exact method. It is a semi-classical approach. It was shown by Lieb and Simon 
that this model will be an exact one for atoms, molecules or solid in general, if the atomic number Z — > oo (4l| (see 
also [12] and [43[ for the very nice review articles). The Thomas-Fermi model in presence of ultra-strong magnetic 
field (when only the zeroth Landau level is populated) in the non-relativistic regime was first used by Kadomtsev 
[44| . It was shown that in presence of strong magnetic field, the electrons move in cylindrical shells with the axis 
directed along the magnetic field. The atoms thus have elongated cylindrical shapes and much more binding energies. 
The Thomas-Fermi model for heavier atoms have also been studied by Mueller et. al., using a variational approach 
[45l | . The relativistic corrections of these calculations was done by Hill et. al. (46j. But in each of these studies an 
extremely strong magnetic field was considered, so that electrons occupy only the zeroth Landau level. 

To develop an exact formalism for relativistic Thomas-Fermi-Dirac model in presence of strong quantizing magnetic 
field at zero temperature (it is assumed that the electron gas within the cell is strongly degenerate) , we assume that the 
magnetic field B is uniform throughout the star and is along ^-direction, i.e., our choice of gauge is = (0, 0, xB, 0). 
Now in the relativistic scenario, the Landau levels of the electrons will be populated if the magnetic field strength B 
exceeds the quantum critical value i?( c )( e ) = m 2 /e w 4.4 x 10 13 G (throughout this article we assume h = c = 1). In 
the relativistic regime the quantum critical value is the typical strength of the magnetic field at which the electron 
cyclotron quantum exceeds the corresponding rest mass energy or equivalently the dc Broglie wave length for electron 
exceeds its Larmor radius. In presence of such strong quantizing magnetic field along Z-axis the electron momenta 
in the orthogonal plane get quantized and are given by p± = (2veB) 1 ^ 2 , where v = 0, 1, 2...., the well known Landau 
quantum numbers. The component along Z-axis varies continuously from —00 to 00, for non-zero temperature, 
whereas, in the zero temperature case, we have the relation: — pp < p z < +pf, where pf is the electron Fermi 
momentum. The phenomenon is known as the Landau quantization. Further, the phase space volume integral in the 
momentum space in this quantized condition is given by 



(2 



I f I f eB v ~ 00 r+°° 

^3 J d 3 pf(p) = j—^ J d Pz d 2 P U(p) = E ( 2 - <W J ^ d Pz f{v,p z ) 



(1) 



The presence of the multiplicative factor 2 — 5 v q is justified by the fact that the zeroth Landau level is singly degenerate, 
whereas all other states are doubly degenerate (which will only be obvious if one solves Dirac equation in presence of 
strong external magnetic field of strength B > B^^). The modified form of spinor solutions of Dirac equation for 
electrons in Dirac-Pauli representation, in presence of strong quantizing magnetic fields are given by 
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where the symbols | and j. indicates up and down spin states respectively, 
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with H v is the well known Hermite polynomial of order v, E v = {p 2 z 



2veB) 1 / 2 , the single particle energy eigen 



value, L y , L z are the length scales along Y and Z directions respectively, e is the magnitude of the charge carried by 
electrons and m is the electron rest mass. 

Now for B > i?( c )( e ), using the inequality p 2 F > 0, the upper limit of Landau quantum number upto which can be 
occupied at zero temperature is given by 
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which is an integer but less than the actual value of the quantity within the third brackets at the right hand side 
and \i e is the electron chemical potential. At zero temperature, the upper limit of v-sum will be ^ max instead of 
oo. The external magnetic field will therefore behave like a classical entity if the strength is less than the quantum 
threshold value and in this region one has to take the standard form of plane wave solution of Dirac spinors. Whereas, 
in presence of ultra-strong magnetic field, the maximum value of the Landau quantum number as mentioned above 
becomes zero. Which indicates that all the electrons occupy the zeroth Landau level have their spins aligned opposite 
to the direction of magnetic field. Following eqn.(l), it is very easy to show that in presence of strong quantizing 
magnetic field the number density for electron is given by 

In the Thomas-Fermi-Dirac model in presence of strong magnetic field in the relativistic region, we replace the 
atoms/ions by Wigner-Seitz cells with varying electron density and assume that A and Z are the mass number and 
atomic number within the cell (A = N + Z, with N, the number of neutrons). To make each cell charge neutral, Z 
must also be the number of electrons within the cell. We further assume that instead of a point object, the radius of 
the nucleus is given by r n = tqA 1 ' 3 with ro = 1.12fm. We will see later that this choice will remove the singularity of 
Thomas- Fermi equation at the origin. The electrostatic potential V(r), felt by electrons satisfy the Poisson's equation, 
given by 

V 2 V{r) = 4nen e (r) - 4iren p (r) (8) 
where n p (r) is the proton density within the nucleus, given by 

%> = S-^-r) (9) 

The second term on the right hand side of eqn.(8) is nuclear contribution. Now the maximum energy of an electron 
at f within the cell is given by 

e v {r) — eV(r) = constant = fi e (10) 
or (pp +m 2 + 2veB) 1/2 - eV{r) = constant = [i e (11) 

In case it is not a constant, the electrons will try to occupy a position within the cell to have minimum energy. This 
will develop an instability in the system. From the above equation we have 

p F = [( Pe + eV(r)) 2 -ml] 1 ' 2 (12) 

where m 2 = in 2 + 2veB. Since we are interested to have electron distribution only outside the nucleus, we discard the 
proton contribution in the Poisson's equation. Further, the potential must satisfy the following boundary conditions: 

rV(r) = Ze for r — > r n (a) 
dV 

—— = for r^r s (b) 

dr 

where r s is the surface value of r (the geometrical structure of each cell is assumed to be spherical in nature). The 
Poisson's equation is then given by 

-^(rV(r)) = — XJ(2 - ^ )[(Me + eV{r)f (13) 
r dr z 7r 
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We now substitute 
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Li e + eV(r) = Ze/-±-^ (c) 
r = [ix (d) 
Then the modified form of the Poisson's equation is given by 

d 2 



= 5>-M(^(*)-tffc 2 ) 1/2 (14) 
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with 



1/2 



and 



m u fi 



From (c) and (d), it is therefore quite obvious, that the radius of each spherical cell, at the crustal region, decreases 
with the strength of magnetic field and the change is oc B~ x l 2 . The squeezing of the Wigner-Seitz cells in presence 
of strong magnetic field is analogous to the well known magneto striction phenomenon observed in classical magneto- 
static problem. Further, the right hand side of the final form of the Poisson's equation (eqn.(14)) must be real. It 
requires, that the inequality <Pqx < 4>(x) must be satisfied. Which after using cj>o and m„, gives 
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This equation indicates that the upper limit of the Landau quantum number of the levels occupied by the Wigner-Seitz 
electrons depends on its position within the cells. Or in other wards, the value of v max at some point x within the cell 
also depends on the strength of electrostatic potential at that point (of course it is assumed that the magnetic field is 
constant throughout the star). Since we have assumed a finite dimension for the nucleus, the problem of singularity 
at the origin will not appear here (origin is actually excluded in our numerical calculation, see also |40|). Since the 
minimum value of v max is zero for ultra-strong magnetic field case, we have a maximum value of atomic radius, given 



by r n 



fIX n 
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< 



2eB 
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The value of x max , obtained from the numerical solution of the Poisson's equation, must satisfy this inequality. 
Further, it is very easy to show from the conditions (a) and (b), using (c) and (d), that the initial and the surface 
condition for V(r) are given by 



(x)\x=x n = 1 and 



d<t>{x) . 
dx 1 



(17) 



respectively. These conditions are identical with the field free case [3l|. While solving the Poisson's equation nu- 
merically using the two boundary conditions (a) and (b) (In the numerical calculation we have followed the standard 
4-point Runge-Kutta technique with shooting method at the surface), we found that it is absolutely necessary to 
incorporate all other extra conditions appearing in this particular case (eqns.(15) and (16)). The surface condition 
must also satisfy 
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r n e (r)dr — Att/i" 



x n e (x)dx 



Since the integration range is from r n to r s (origin is avoided), the serious problem of singularity associated with 
Thomas-Fermi-Dirac equation, at the centre of the cell will no longer be there. Therefore, in this particular case, it 
is also not necessary to follow the numerical methods prescribed by Feynman, Metropolis and Teller (47| . 

The corresponding expression for the Poisson's equation for B = in the non-relativistic regime is given by [48j 
(see also [3l[ and (49j 1 
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dx 2 



^3/2 

7m 



where 



0) = 1 and 4>'{x s ) 



<j)(x s ) 



and 



9tt 2 \ 
128Z J 



1/3 



ao 



with ao = Ti 2 /{me 2 ), the Bohr radius. This equation has a singularity at the origin. 

To obtain the variation of <j>{x) with x within a typical Wigner-Seitz cell, we have solved the Poisson's equation 
(eqn.(14)) numerically considering all the necessary conditions as mentioned above along with the boundary conditions 
(a) and (b). In fig.(l) we have plotted <p(r) as a function of r (in A), the radius of the cell, for three different initial 
values for </>' and magnetic field strength B — 10 14 G. In fig. (2) we have plotted the same quantity but for three 
different magnetic field strengths. It is quite obvious from the curves of fig- (2) that the surface of the cell is reached 
later for low magnetic field strength compared to stronger values. This is consistent with the well known magneto 
striction phenomenon. 



6 



3. EQUATION OF STATE 



To obtain the equation of state of such low density crustal matter in presence of strong quantizing magnetic field 
in the relativistic region, we first evaluate the kinetic pressure, given by 
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This momentum integral can very easily be evaluate and is given by 

'PF + ipF+ml) 1 ' 2 



2tt 2 



v=0 



p F (p 2 F + ml) 1 ' 2 - m 2 v \n 



(18) 



(19) 



where pf = Pf(x s ). The above expression, therefore gives the cell averaged kinetic pressure. To obtain the numerical 
values, it is necessary to express pf(x s ) as a function of (f)(x s ) using eqn.(21) as given below (cither the correct 
expression, particularly for numerical evaluation or the approximate one for the analytical result). For the exact 
expression, as given in eqn.(19), the functional relation is extremely complicated and has to be obtained numerically. 
To get a simple mathematical expression for kinetic pressure, we consider the ultra-relativistic case, which gives 



4tt 2 



(20) 



u=0 



Now to express the Fermi momentum as a function of (f>(x s ) in the ultra-relativistic limit, we consider 



PF = [(^ + eV(r)) 2 -ml}^ 



2„4 



fix s 



1/2 



which is the correct expression. On putting this correct expression for Fermi momentum in eqn.(19), the exact 
functional dependence of kinetic pressure on 4>(x s ) can be obtained. The approximate form may be obtained by 
neglecting to„, and we have 



p F « Me + eV(r) 



Substituting the conditions (c) and (d), we have 
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p F « Ze 



eB 
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In this approximation, the electron number density is given by 

j-) Z-'max 

He- 



v=0 



■, <KXs) 
LLX S 



(21) 



(22) 



(23) 



Whereas, the usual form of mass density or rest energy density of the crustal region of neutron stars, contributed by 
the massive nuclear part, is given by 



p(x s ) = e(x s ) 



3AniB 



(24) 



Here, the nuclei are assumed to be static in nature and are not affected by strong magnetic field. Whereas, electrons 
within the cells are considered to be almost free and dominates the kinetic pressure part. On eliminating x s from 
eqn.(24) and either from eqn.(19) or eqn.(22), we get the equation of state of low density crustal matter in exact or 
approximate form respectively. 
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4. KINETIC ENERGY FOR ELECTRON GAS 



In this section we shall evaluate the energy contribution from the electron part within the cell. Although it is 
small enough compared to the rest energy from nuclear part in the non-relativistic region, we have noticed that it 
has a significant amount of contribution in the relativistic scenario. The energy contribution from the electronic 
part consists of four main parts, they are: (i) the kinetic energy part, (ii) the electron-nucleus interaction part, (iii) 
the electron-electron direct interaction part and (iv) the electron-electron exchange interaction part. There can be 
another important contribution, the correlation part. Now the correlation energy is not a quantity with physical 
significance. It actually gives the error incurred in making a fairly crude approximation. In the Hartree zeroth order 
approximation, the TV-electron distribution function factors into a product of N one-electron distributions. In the 
Hartree-Fock wave function (Slatter determinant form) correlations are introduced in the first order approximation. 
In the evaluation of actual correlation term, higher order approximations are incorporated and the exchange term is 
excluded. Since the exchange term is absorbed in the many body correlation terms, in the density functional theory, it 
is conventionally known as the exchange-correlation term and obtained with some approximation of electron density. 
In this article, however, we have not taken the correlation part into account. In some future communication we shall 
bring this extra term in the energy expression and evaluate numerically the equation of state of low density crustal 
matter. 

Now the contribution from electrons to the kinetic energy part is given by 



fa p fj fPF 
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Evaluating the integral over p z , the exact expression for kinetic energy is given by 
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where the exact form of pf is given by 



p F (x) = 
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(27) 



Again this will give an extremely complicated functional dependence on <p{x). In addition, since u max {x) depends 
on the radial coordinate, the evaluation of Eke by analytical integration is just impossible. Therefore, numerical 
technique has to be followed to evaluate the integrals along with the sum over v. Numerically fitted functional forms 
for (j){x) and the corresponding v max {<i){x)) are used to evaluate the kinetic energy numerically. To get a simpler 
expression, we go to the ultra- relativistic limit, then 
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To obtain an approximate result, we use the relations 
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Then integrating by parts the Poisson's equation, considering the surface condition and finally with the trivial relation 
(we further assume that ^ max is a function of x s , instead of a function of x) 



^(2- 5 v0 )v = v m „(i/ mm + 1), 



v=0 
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we have 



E 



KE 



Z 2 e 2 
2fi 



f>'(x n ) } + (X 3 S ~ ^)^ max (i/ max + 1) 



For ultra-strong magnetic field, i/ max = 0, then we have 

&KE 



2(1 
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On the other hand, from eqn.(27) if we consider 
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(31) 
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then in the ultra-relativistic approximation we have the same result as shown in eqn.(32). The cell averaged electron 
kinetic energy density is then given by 



Eke %Eke 3Eke 



The corresponding expression for B = in the non-relativistic regime is given by [48 
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5. INTERACTION ENERGY:ELECTRON-NUCLEUS 



Next we consider the three possible types of interaction term. Let us first consider the electron-nucleus interaction 
term, given by 



E„ 



-Ze 1 



3„ "e 



A-KZe 2 (i 2 I xdxn e {x) 



(35) 



where the exact expression for electron density n e (x) is given by eqn.(7) and the Fermi momentum pf{x) is given by 
eqn.(27). Then we have after expressing pf(x) as a function of <)>{x) 



E en = 5(2 - 5 u0 )dx(4>(x) 2 - cgx 2 ) 1 ! 2 



(36) 



This integral can not be evaluated analytically. To obtain the exact values for E en , as before, we need <fi(x) as a 
function of x and ^ max as a function of <f)(x) within the integration range. To get a simple form, we consider the 
ultra-relativistic limit. Then we have 



2 C 2 r x 



E„ 



Z z e 



y^(2— 5 v0 )dx(f>(x) 



Using the approximate form of the Poisson's equation, given by 

"man 



we have (again assuming v max (x s ) instead of a function of x) 

Z 2 e 2 r x s 
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(39) 
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We finally have, after integrating by parts and using the surface condition for 4>(x) 



Z 2 e 2 



4>(x s ) 



<fi(x n ) 



The sum of electron kinetic energy and the electron-nucleus interaction energy is then given by 



Eke + E en — — - 



Z 2 e 2 



2(i. 



4>(x s ) 



and the cell averaged electron-nucleus interaction energy density is given by 

3E en 3E 

en. 

V 4nr 3 4tt^ 3 x 3 
The corresponding zero field value in the non-relativistic regime is given by [46 

Z 2 £ 2 rx a 



(40) 



(41) 



(42) 



E — — - 



^ 2 (x)x^ 2 dx 



V Jo 

M Jo 

Z 2 e 2 (<j>{x s ) 



0'(O) 



6. INTERACTION ENERGY:ELECTRON-ELECTRON DIRECT INTERACTION 

Next we consider the electron-electron interaction. Let us first evaluate the direct term. It is given by 

Eg = \e 2 J d\n e {r) J d 3 r'n e (r')-^j (43) 

Assuming r as the principal axis and 9 is the angle between r and r', we have d 3 r — 47rr 2 dr, d 3 r' — 2irr' 2 dr' sin9d8 
(we have assumed that the vectors f and r' are on the same plane) and \r — r'\ = (r 2 + r' 2 — 2rr' cos 6 1 ) 1 / 2 . The limits 
for both r and r' are from r„ to r s and the range of 9 is from to n. Let us first evaluate the angular integral, given 
by 

I(r,r') 



lo (r 2 +r' 2 - 2rr'cos6») 1 /2 
It is straight forward to show that I(r, r') = (r + r') — \r — r'\. Then we have 

£W = 47r 2 e 2 / rdrn e {r) I r'dr'n e (r')[(r + r') - \r - r'\ 



(44) 



Now from this equation it is trivial to show that the quantity within the third bracket will be 2r' for r' < r and for 
the opposite case it will be 2r. Then the above expression, for direct interaction, reduces to the following simple form: 



Eg = 8e 2 7r 2 | J rdrn e (r) J r' 2 dr' n e (r') 
+ I r 2 drn e (r) I r'dr'n e (r')\ 



(45) 



To obtain Eg, we have evaluated separately the integrals appearing in the above expression. They are given by 

Ii = rdrn e (r) x l2{f) = [J 2 / xdxn e (x) x l2(x) 
h = I r 2 drn e (r) x 7 4 (r) = fi 3 / x 2 dxn e {x) x h(x) 
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where 



pT f-X 

h{r) = / r' 2 dr'n e (r') = M 3 / x' 2 dx'n e (x') = h{x) 
h(r) = I r'dr'n e (r') = f x'dx'n e (x') = h{x) 

J r J x 



It is easy to verify that the only integral which can be evaluated analytically is I±{x). It is therefore necessary to 
obtain I\, I 2 and I 3 numerically. Therefore, the exact analytical express for Eee can not be obtained. However, we 
can have some approximate analytical results. As for example, 



I 2 {x) = /i : 



rx ^ 

/ x' dx'n e (x') 

J X n 



e 3 ZBfi 
2tt 2 

e 3 ZBfi 
2ir 2 

e 3 ZBii 2 
2ir 2 



2 rx 



^(2-^ o )(0 2 -^' 2 ) 1 /Vcix' 



. . .g^ ^ 6 v0 )<f>(x')x'dx' 

r ^(x'ydx' 

J Xn 



47T 



[{x(f>'(x) - X n (j)'{x n )) - (</>(x) - <f)(x n ))] 



Similarly, we have approximately 

h « M 2 

Z 2 e 2 eB f ia 



xdxn e (x)— [(x(j)'(x) - x n (j>'{x n )) - (4>{x) - </>(x n ))] 

47T 



(j)(x)dx[(x(j)' (x) - x n <p'{x n )) - (<j>(x) - (f>(x n ))} 



To evaluate this integral, we split it into four integrals, given by 

= f ' (f>(x)(f>' (x)xdx 

i[ 2) = = r <p(x)dx 



i[ 3) = HimTdx 

J x n 



The approximate results for l[ 2 ^ (=l[ 4 ^) is straight forward and is given by 



Similarly, 



l[ 2) = l[ 4) = cj>\x s ) - 4>'{x n ) 



= \ [x s <j){x s ) 2 - x n (f)(x n ) 2 ] i and l[ 3) = 



Then we have approximately (again assuming that f max is a function of x s ) 



Z 2 e z neB 
8tt 3 



- {x s cj)(x s ) 2 - x n (f>(x n ) 2 } 



- 4>'{x n )x n {<t>'(x s ) - 4>'{x n )) + 4>(x n ) (4>'(x s ) - 4>'{x n )) 



(46) 



(47) 



(48) 



On putting this result in eqn.(45), the approximate form of direct part of electron-electron interaction energy contri- 
bution from the integral I\ is given by 



2/j. 



{x s (j)(x s ) 2 - x n (j)(x n ) 2 } 



4>'{x n )x n {4>'(x s ) - 4>'{x n )) + 4>{x n ) {4>'{x s ) - 4>'(x n )) 



(49) 
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where a — v mH (i/ mm + 1). Similarly the integral I4 is approximately given by 

eBZe 2 /i [ x °" 



2ir 2 



[ ' ]T (2 - 6„ )<t>(x')dx' 



eBZe 2 /i f Xs „ , , 



2tt 2 

eBZe 2 fj, [ x ° fAl 
2tt 2 



4>"{x')dx' 

rv(* s w(x)) 

J X 



Similarly the integral J 3 may also be approximated by the integral 

_ fieBZ 2 e 2 <- x " 



8tt 3 



This integral may be broken into two parts 



2( 2 - 5 M )^{x)xdx{4>'{x s ) - 4>'{x)) 



»«ai /.X° "max 

^ (2 — <5 1/0 ) ( / ) (a ; )a ; rfa ; an d = / ^ (2 — 8 v q)4>' {x)4>{x)xdx 



v=0 



r(l) 
i 3 



The first one is approximately given by 

w / (f>"(x)xdx 

J x n 

= {x s (j)'(x s ) - x n (j)'(x n )) - {4>{x s ) - 4>{x n )) 
and the approximate value for the second one is 

4 2) ~ / 4>'(x)(f>"(x)xdx 

= x s (4>'(x s )) 2 - x n (4>' (x n )) 2 
Combining these two results, we have the approximate form of I3 as given below: 



/ 3 



/ieBZ 2 e 2 
8tt 3 



{(x s (j)'(x s ) - X n 4>'{x n )) ~ {4>{X S ) - 4>{x n ))} - 

{x s {<j>\x s )f -x n (<j>'(x n )) 2 } 



(50) 



(51) 



(52) 



(53) 



(54) 



(55) 



The approximate form of direct part of electron-electron interaction energy, coming from the integral I3 may then be 
obtained by substituting eqn.(55) into eqn.(45) and is given by 



Z 2 e 2 

2/2 



(56) 



{(x s <t)'(x s ) - X n (j)'{x n )) - (<p(x s ) ~ 4>{x n ))} - 

{ x s (^(x s )) 2 -x n (^(x n )) 2 } 

Finally adding eqns.(49) and (56), we have the approximate form of direct part of electron-electron interaction energy 

Eg = + E(?M (57) 



The cell averaged direct interaction energy density is then given by 



Ad) = 

c ee 



3E, 



(d) 



V 47rr 3 4irn 3 x 3 



(58) 
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The corresponding direct part of e-e interaction term for the zero field case in the non-relativistic regime is given by 
Eg = — r x 1 / 2 ^/ 2 {x)dx{- I* x' 1/2 <j?l 2 {x')dx' + f S x'- 1,2 ^/ 2 {x')dx' 

2 M io Uio Jx 

It is straight forward to evaluate these integrals using the Poisson's equation and the surface condition. Then we have 

J*yV /2 (^)-^'(o)" 

The exact values for EeV can only be obtained numerically, knowing (j>(x) as a function of x within the range x n to 
x a . This is true even for the field free case. 



E W = I™ 
2 [i 



7. INTERACTION ENERGY:ELECTRON-ELECTRON EXCHANGE TERM 

Next we shall consider the exchange term. The exchange energy integral corresponding to the ith. electron in the 
cell is given by 

E( eT ] = T E / d\d\'-^M^ 3 (?)^Mr') (59) 
1 j J \r-r\ 

where the spinor wave function is given by eqns.(2)-(4) and 'tp(r) = ip'{r)"fo, the adjoint of the spinor and 70 is the 
zeroth part of the Dirac gamma matrices 7^. Now it is very easy to show that for t = t' 

— —* 2771 

^i(^)ipi(r') = exp[-i{p y (y - y') +p z (z - z')}] 

)I^ hPy (x')} (60) 



Similarly, we have 



— — * 2 777, 

^j(r')^j(r) = T T exp[i{p' y (y-y')+p' z (z- z')}] 



{I»';p' v (x)Iv'; Py {x') + IW-l;p' y (x)Iu'-l;p' y (x')} (61) 

When these two terms are combined, we have after replacing the sum over j by the integrals 

/+oo r+PF 
d P'y / d Pz 
-oc J -~pf 



exp[-i{{p y - p' y ){y - y 1 ) + (p z - p' z ){z - z')}] 

[{Iv,Py {x)Iu;p y (x') + L-l; Py (x)I„-i ;Py (x' )} 

{I v ,, ply {x)I vl , K {x') +I v ^ v>Py (x)I u ,- Uply (x')}} (62) 
It is possible to evaluate the integrals over y' and z', given by (50j 

j dy'dz'— — -■ exp[-i{(py-p y >)(y-y') + (p z -p z >)(z-z')}} 

-00 J —00 \T T I 

A-7T 

= —expi-Klx-x']) (63) 
where K = [(p y — p y f ) 2 + (p z — Pz') 2 ] 1 ^ 2 ■ Then the integral over y and z is given by 

4-oo p-\-oo 

/ dydz = L y L z 

oc J — oc 
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Since the exchange integral is multi-dimensional in nature with an extremely complicated form of integrand, it is just 
impossible to evaluate the exchange energy analytically. However, a simple form can be achieved if all the electrons 
occupy only the zeroth Landau level (i.e., i/ max = 0). This special case have been discussed later in this article. In 
the numerical evaluation of the above exchange energy integral, we have put p z — pf, the electron Fermi energy and 
for the sake of convenience we make the following substitution: 



eB 



P y =Py- Py and P z = p' z - p z 



Since — oo < x < +00, —00 < x' < +00 and —00 < p' y < +00, the limits of the new variables will also remain same. 
Further, we have used — pf < p' z < +Pf and as consequence, we have < P z < 2pp- Although the ranges of the 
variables x, y, x' and y' are from — r s to — r„ and then from +r„ to +r s , for the sake of convenience we have considered 
the limits from —00 to +00, assuming that the electron density vanishes inside the nucleus and also outside the cell 
surface. 

The exact values for the exchange integrals are evaluated numerically using Monte-Carlo technique to obtain 
electron-electron exchange interaction part. In fig. (3) we have plotted the exchange energy (in MeV) for various 
values of electron Fermi momentum (in MeV) . We have fitted numerically the exchange energy as a function of Fermi 
momentum and is given by 



E ee (ex) = exp(ap F 



where the parameters E^ ex> = 0.043 and 0.062 in MeV and a = 0.409 and 0.42 in MeV" 1 for B = 10 14 G and 10 16 G 
respectively. It has been observed that the minimum value of pf for which E^te^ is non-zero increases with the increase 
in B. The qualitative nature of the curves are exactly identical. However, E^t^ increases with B for a given pp. 

We shall now consider the Thomas-Fermi-Dirac model for v = 0. We assume here that the strength of magnetic 
field is extremely strong, so that electrons occupy only the zeroth Landau level {v maK = 0). In this scenario, the 
mathematical derivations are much more easier than f max 7^ case. In this approximation 



e F = (p 2 f + m 2 ) 1/2 and 
Further, we define the electron chemical potential 



eB 
2^r 



-VF 



Which gives 



He = e F {r) - eV(r) 



eB 



1/2 



— eV(r) — m = constant 



— {2m(eV(r)+n e )} 



1/2 



(eF(r)+ Me ) 



2m 



1/2 



(64) 



(65) 



(66) 



Substituting this expression in the Poisson's equation (eqn.(8)) and discarding the nuclear contribution, we have 



p B 

y 2 V = ATTe^r\2m(eV(r) + /i e )] 1/2 
Following the same procedure as we did for v ^ 0, we get 



(eV(r) + /i e ) 



2m 



1/2 



dx 2 



1/2 



1 + A- 



where 



A = 



Ze 2 
2m/i ' 



r = fix where [i — 



1/2 



^1/5^2/5 



(67) 



(68) 



2 3/5 e 4/5 B 2/5 TO l/5 



The initial and the surface conditions remain same as we have for t^ max 7^ case. Further, the form of this equation 
is such that there is no singularity at the origin. Therefore, in this case we need not have to integrate from r n (we 
can assume safely that the nucleus is a point object), or do not have to use the numerical prescription as discussed in 
reference 47 1. 
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Almost the same algebraic procedure is followed to obtain Thomas-Fermi-Dirac equation for the relativistic case 
with B = [H. In this case 



e = [p 2 + to 2 ] and n e = —^p% 



Here, we have 



Which gives 



fi e = e F (r) - eV(r) - m = [(37r 2 n e ) 2/3 + to 2 ] 1/2 - eV(r) - m = constant 



2m(p e + eV(r)) 
3^ 



3/2 



(eV(r)+[i e ) 



2m 



3/2 



Then substitution of these expressions in the Poisson's equation, gives 



dx 2 



63/2 

B i/2 



3/2 



nuclear contribution 



(69) 



(70) 



(71) 



where 



3tt 



1/2 



Z cr = ( ^3 ) and /i = 



(37T) 



2/3 



s 2 2 7/3^1/3 



The initial and surface conditions are again same as for B ^ and ^ max ^ 0. Unfortunately, it has singular nature 
at the origin. Which is removed artificially by taking the lower limit of r or ^-integration as r„ or x n respectively or 
following reference [47| (see [I(|). Similarly for the non-relativistic regime with B ^ but i/ max = 0, it is very easy to 
formulate all the above equations. Here 

( \ P%(r) i eB 



Further, 



It gives 



2to 



eV(r) = constant 



B 



n e = — pm^ + eVir))} 1 / 2 
2ir z 

Which after substituting in the Poisson's equation, we have 

d l± = ^1/2^1/2 

dx z 

This equation again has no singularity at the origin. Further, in this case 

^2/5^1/5 



(72) 



2 2/5 e l/5 B 2/5( 2TO )l/5 



but the boundary conditions again remain as usual. 

Let us now consider the non-relativistic case with B ^ and t-Wx ^ 0. In this case the energy eigen value is given 
by 



2m \ 2 1m 



Putting p z = pf, we have 



P F 



2 to 



1\ eB 



[i e = — — + I v + - ] eV(r) = constant 



2/ to 
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Hence, 



eB 
2ir 2 

eB 



E 



Pf 



= ^Ei 2TO (Me+^W)-(2^+l) e S} 1 



/2 



Following the same procedure as we did for the relativistic case with J/ max / 0, the Poisson's equation reduces to 

dx 2 



,1/2 



1/2 



(73) 



where 



^1/5^2/5 



2 2/5 e 2/5( eB )2/5( 2m )l/5 



and 



veBp 
mZe 2 



In this particular scenario, by inspection one can realize that the necessary condition to be satisfied to have physically 
acceptable solutions is that the quantity within the square root on the right hand side of eqn.(73) should not be 
negative, i.e., 

mZe 2 (j>(x) 1 
~ 2 



x<t> v > which gives f ma x < 



eB/ix 

The above form of the Poisson's equation with both B and i/ max ^ 0, has no singularity at the origin. The boundary 
conditions are same as for the relativistic case. The above condition imposed on the upper limit of v has to be checked 
at every steps of numerical integration. Since f ma x > 0, we have 

mZe 2 4>(x) ^ 1 



eB/ix 



The wave function ip{r) cx I v - tP (x) for the electrons with 7„ ;p , given by eqn.(5), is quite complicated. The primary 
reason is the presence of Hermite polynomial of non-zero order. The exchange energy in this case also can not be 
evaluated analytically. The expression is as complicated as we have derived for the relativistic case. 

On the other hand, with zero magnetic field for non-relativistic scenario, the exchange energy is given by [5l| (see 
also jH) 



(p 



■In 



Pf + Pz 



Pf - Pz 



2pi 



(74) 



Substituting p z — pp, we have 



Pf 



(75) 



For the sake of completeness we shall now briefly discuss the evaluation of exchange energy for both the non- 
relativistic and relativistic cases in presence of ultra-strong magnetic field, i.e., B ^ and ^ max = 0. For both these 
cases, semi-analytic expressions can be obtained. The exchange energy in the non-relativistic regime is given by 



3 = 1 J 



•3„j3 r /_ 



1 



-ip* (r)tpj {rjip* {r)ipi (f 1 



For the zeroth Landau level, since Hq(x) — 1, the wave function ip(r) for electron is given by 



1 



(L v L z y/ 2 V n J 



eB\ 



1/4 



exp 



eB 
~2~ 



Py_\ 
eB) 



exp[i(p y y +p z z)] 



(76) 



(77) 



Here also we can replace the sum over j by the integrals J J L y L z dp' y dp' z . Writing d 3 r' — dx'dy'dz' . and following 
Lee 1501, we have 



/ 



1 



dy'dz'—. — — exp[~i(p y - p' )(y - y 1 ) - i(p z - p' z )(z - z')] 



47T 

2K 



exp {—K x — x' |) 



(78) 
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where A = \(p y - p') 2 + (p z - p' z ) 2 }^ 2 . 



Similarly for d 3 r = dxdydz, the integral J dydz = L y L z . Then we have 

Et x) - \ (^) 47re 2 J dp' y dp' z dxdx' ± exp (-A | x - x' 



exp 



■t{(*-%)'+(*-£)'+(*-S)'+(*-£)'1 



(79) 



To evaluate the integrals over x and x', we change the integration variables to X and Y, where X = x — x' and 
Y = {x + x')/2. 
Now 



J — ( 



rfXexp(- J ft' | X |) exp 



-X' 



2vr 



1/2 



exp 



if 2 \ 



2eB , 



erfc 



A' 



(2eB) 1 /2 



(80) 



where erfc(x) is the complimentary error function. 
Then we have 



Ei e e x) = e 2 B J ±dp' y dp' z dY 



2tt 

~&~b 



1/2 



exp 



( K 2 



\2eB 



erfc 



(2eB) 1 /2 



exp 



eP 



The y integral is given by 



f 

j — i 



dYexp 



4F 2 + 



£^ ( 2 y Py+Py 



2 P 2 y ^ 2p' 2 4 Py Y Ap' y Y 



+ 



e 2 B 2 e 2 B 2 eB 



eB 



eB 



2eB 



1/2 



(81) 



(82) 



Then we have after changing the integration variables from p' y and p' z to P y = p y — p' y and P z = Pz — p' z 



E^ = e 



! 7T j dPydP Z 



(P 2 + P|)V2 exp \2eB 



( P 2 

I ± Z 



rfc 



■p2 + p2\ 1 / 2 
± y * ± z \ 



2eB 



(83) 



where the limit of P y is from — oo to oo and P z is from to 2pp for p z = pp. 
Again putting P y = P z tan#, we have 



e 7r 



2, n . 



dP z 



it/2 



sec d9 erfc 



Pz 



(2eP)V2 



sectf exp 



2eB 



(84) 



This is the form of semi-analytic expression for the exchange energy with f max = in the non-relativistic regime. Fur- 
ther simplification of this expression is not possible. Therefore, these double integrals have been evaluated numerically 
as a function of Fermi momentum pp. The fitted functional form of E^T^ is given by 



Et x) = a[l - exp(-0p F )] 

where the parameters a and (3 vary with magnetic field strength B and are shown in the following table. 
Table-I 



(85) 
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B (Gauss) 


10 14 


10 15 


10 17 


a (MeV) 


0.568 


1.796 


17.909 


(3 MeV" 1 


3.412 


1.067 


0.109 


7 


0.506 


0.527 


0.658 


C 


0.973 


0.870 


0.386 


X s 


3.096 


3.170 


4.404 


Ts (A) 


0.402 


0.203 


0.123 


V s 


-0.938556 


-0.937365 


-0.936123 


<t>0 


1.633 


1.651 


1.944 


e 


2.097 


2.071 


1.755 




0.213 


0.204 


0.031 


P (gm/cc) 


72.79 


572.29 


962.14 



Unlike the relativistic case here one can see from the fitted functional form, that the exchange energy saturates to 
some constant value a(B). Now, if we include the exchange part separately, then in Thomas- Fermi-Dirac model the 
electron Fermi energy is given by, 



u, = — ecj) — E^l x \p F ) — constant 

2m 

Rearranging the above equation in the form (see also [52|). 

-Ppf 



2m 



ae 



(86) 



(87) 



where fi* = fi + a is the modified form of Fermi energy of the electron, one can express Fermi momentum as a 
function of //* + e<f>. The numerically fitted functional form is given by a simple power law, 



PF = C(fj,* + 



(88) 



where, C and 7 are constant parameters for a given magnetic field strength. In Table-I above, we have shown the 
variation of C and 7 with the magnetic field strength B. The variation of <j)(x) with x for a given magnetic field 
strength is given by the numerically fitted functional form (the solution of the Poisson's equation is fitted numerically) 



<j>(x) 



1 + exp{£(a; - xo)} 



(89) 



where, 0o, £, %o are constant parameters for a given magnetic field strength. The variation of these parameters with 
magnetic field strength are also shown in Table-I. In presence of strong quantizing magnetic field, the variation of <fi 
with x is entirely different from the non-magnetic case. The variation is more or less like the radial distribution of 
matter in neutron stars. Further, we use 



dx 



(j)(x s 



We have also shown the variation of v s with B in Table-I. 

We next consider the relativistic case with z/ max = 0. Since for v = 0, we have H u (x) — 1 and H v -\{x) — 0, the 
appropriate form of exchange energy is given by 



E. 



{ex) = e 2 J ... J ' dp' y dp' z dxdx'^~^{-K\x-x'\) 

h;p y (x)I ;p v {x')Io-p' (x)I ;p' (x') 



(90) 



where 



eB 

7T 

X = 



1/4 



exp 
eB) 



-eBX- 
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Py = Py- Py 
Pz = P' z - Pz 
Pz = PF 

E = (p' z 2 + my 

We again get a similar type of semi-analytic expressions as shown in eqn.(84). The physical quantities in this scenario 
are when evaluated numerically, can be fitted by the same type of functional forms, as shown by eqns.(85) and (87)- 
(89). The numerical values for the parameters are more or less same as shown in Table-I. The qualitative nature of 
the dependence of the parameters on the magnetic field remain almost unchanged. A complete numerical analysis of 
this formalism, along with the numerical estimate of equation of state of the crustal matter of strongly magnetized 
neutron stars will be presented in a future communication. In that correspondence we shall also make comparative 
studies of various models and approximations. In this article we have presented very briefly some of the numerical 
estimates to give a feeling of our formalism. 



8. THOMAS-FERMI INDUCED CHARGE DENSITY 

Finally we shall discuss the appearance of Thomas-Fermi induced charge density inside the cells. The total charge 
density within the system is given by 

P{r) = Pext(r) + Pind(r) (91) 

where p- ln d{r) is the induced charge density and is given by the fundamental equation of non- linear Thomas- Fermi 
theory: 

And = -e[n e (fj, e + eV(r)) - n e (/x e )] 
= -e- — eV{r) 

dfl e \V(r)=0 

gE( 2 -M ( ^_% )1/2 ^) 



= — e 

>7T- 

v=0 



= E ( 2 - S*M + mlf' 2 eV(r) (92) 

^ PF „ 
v— 

In the present case, we assume that the electric field V(r) is a slowly varying function of coordinate r. In actual 
practise, to obtain the charge density, one has to solve self-consistently the Dirac equation in presence of an external 
magnetic field (in this case the strength B > £?( c )( e )) and the electrostatic field V(r) and the Poisson's equation 
satisfied by V(r). In the exact scenario, the definition of electron density is given by 

n e (r) = ^(r)i>(r) (93) 

Hence one can obtain the total charge density and V(r). However, the method is extremely complicated, even 
numerically. In Thomas-Fermi approach, we actually do not solve the Dirac equation, instead, assume that V(r) is 
changing slowly with r and get an approximate result. 

Now taking the Fourier transform of the last relation as given in eqn.(92), we get 

Pind(g) = x(q)V(q) (94) 
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Then the Thomas-Fermi dielectric constant is given by 



47T 

= 1 - —x(q) 
q 



i 



(95) 



It is obvious that x(q) (and hence e(q)) is independent of q. Now for the non-relativistic case (m u 3> Pf) 
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On the other hand, in the relativistic scenario (jpp 3> m u ), we have 
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Which gives 
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It can be shown that in presence of ultra-strong magnetic field for both the non-relativistic and the relativistic cases 

(^max=0), 



I?* 

" n 



el ( , P% 
TO+ — 

Pf \ 2to 



To obtain the screened coulomb potential, we use the well known relation 

1 



V(q) = ^V^(q) 



where 



Hence, the screened coulomb potential is given by 



r 



Q 



V ind (r) = — exp(-k r) 
r 



(99) 



(100) 



(101) 



(102) 



Since fc is a function of magnetic field strength, the above equation gives the screened coulomb potential in presence 
strong quantizing magnetic field. In the numerical evaluation of the screening length in presence of strong magnetic 
field, one can use any one of these expressions as given above. 

In this context, we must mention that in a very recent work, Shabad and Usov have investigated the effect of strong 
magnetic field on coulomb potential [53]. It has been shown that the coulomb potential gets modified significantly 
in presence of strong quantizing magnetic field. This is, unlike the Thomas-Fermi model, is an exact field theoretic 
approach. In this work the modified form of vacuum polarization (grossly speaking, this will give a modification of 
the screening length) in presence of strong quantizing magnetic field has been considered. We expect that a lot of new 
results can be obtained if one incorporates these results in Thomas- Fermi model calculation for the crustal matter of 
magnetars. In particular the interaction terms (electron-nucleus, electron-electron direct and exchange terms) will be 
modified significantly and thereby affects the equation of states of this low density matter. 
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9. CONCLUSIONS 

In this article we have developed the formalism for relativistic version of Thomas-Fermi-Dirac model in presence 
of strong quantizing magnetic field. The formalism is applicable to the outer crust of magnetars and also to strongly 
magnetized white dwarfs. 

We have compared our results with several other cases, e.g., the well known non-relativistic model with zero 
magnetic field, field free relativistic case, non-relativistic model in presence of strong quantizing magnetic field for 
both v max / and u max = 0. 

We have noticed that in this formalism, to solve the Poisson equation numerically it is necessary to include a few 
more conditions, which were absent in the usual field free non-relativistic model or in presence of ultra-strong magnetic 
field (V max = 0). 

To remove singularity at the origin, we suggest, following [40j, to use finite dimension for the nuclei. It has also 
been noticed that unlike other scenario, one extra condition appears in the non-relativistic regime with B / and 

fmax ^ 0. 

We have also given an approximate method to get an estimate of the induced charge within each cell and thereby 
obtain the variation of screening length with magnetic field strength. 

In our model, the Wigner-Seitz cells are assumed to be spherical in nature and found that the radius of each cell 
decreases with the increase of magnetic field strength. The variation is given by ~ B^ 1 / 2 . 

The formalism is of course not applicable to the inner crust region, where the matter density is close to the neutron 
drip point, some of the neutrons may come out from the cells. In some future communication we shall present a 
modified version of this formalism appropriate for the inner crust region. 

We have assumed that all the electrons within the cells are moving freely, i.e., they are not bound in any one of 
the atomic orbitals. In reality, it may happen that the electrons at the vicinity of the nucleus in a cell have negative 
energy. These electrons, therefore can not be treated as free. It is therefore absolutely necessary to get the total 
energy of an electron as a function of its position (r or x) within the cell from the numerical solution of the Poisson's 
equation and the expressions for kinetic and various form of interaction energies. We expect that very close to the 
nucleus, the electron energy will be negative and for a particular value of x (= r//i) (which may be a function of B) 
it will become zero (quasi-free electrons) and then becomes positive. If it is found so, then we can not assume that all 
the Z-electrons in the cell are participating in statistical processes. On the other hand, if we consider the expression 
for electron energy as given in eqn.(86), then from the physics point of view all the electrons will become free (energy 
is always positive). Whereas, if we consider 

/i = kinetic energy — e<fi = constant, 

then we may have bound, quasi-free and free electrons within the cells. The presence of free electrons in the compressed 
cells in a dense medium is popularly known as statistical ionization (see reference [HJ] for a detailed discussion). 

To conclude our results, in the following we have given in tabular form the variation of Fermi momentum, Pressure, 
and various kinds of energy (except the exchange energy part) for electron gas within a typical Wigner-Seitz cell, with 
the strength of magnetic field. 
Table-II 



B/B c 


PF (x s )(MeV) 


P{x s )(McV 4 ) 


E KE {x s ){MeV) 


E en (x s )(MeV) 


E^{x s )(MeV) 


10 5 


2.57 


8.64 x 10 4 


15.95 


46.71 


8.12 


5 x 10 4 


2.58 


7153.02 


18.13 


27.15 


12.15 


10 4 


2.66 


4323.14 


19.58 


4.88 


19.80 


5 x 10 3 


2.76 


681.79 


24.53 


-0.37 


23.82 


10 3 


3.56 


319.78 


27.30 


-7.89 


36.18 


5 x 10 2 


4.82 


62.99 


37.67 


-9.69 


44.55 


10 2 


6.95 


38.06 


44.81 


-12.09 


78.04 


50 


9.13 


31.49 


195.02 


-13.13 


93.95 


10 


13.61 


10.06 


781.47 


-13.85 


104.24 


1 


19.94 


6.28 


1.13 x 10 5 


-35.94 


130.67 



From the above tabular form of data one can see that the electron Fermi momentum and the corresponding kinetic 
energy decreases with the strength of magnetic field. Since the exchange energy has to be subtracted and its magnitude 
increases with the magnetic field strength, we may conclude that the system becomes more and more stable (total 
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energy decreases) with the increase in magnetic field strength. 
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FIG. 3: The variation of exchange energy with Fermi momentum of the electrons. The magnetic field strength B — 10 14 G 
(solid curve) and 10 16 G (dashed curve) 



